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I.  INTRODUCTION 


Measurements  of  rate  coefficients  for  energy  transfer  and  deactivation 

processes  in  HF(Vj)  + DFfv^  = 0)  collisions  have  not  been  reported  in  the 

literature  for  v.  > 1.  The  laser-excited  fluorescence  method  has  been 
1-6  * 

used  to  measure  the  temperature-dependent  quenching  rate  coefficient 

for  (v— v + v — R)  processes  in  HF(Vj  = 1)  + DFfv^  = 0)  collisions.  These 

measurements  indicate  rapid  coupling  between  the  vibrational  and  rotational 

14  6 

motions.  The  laser-excited  fluorescence  technique  has  been  used  ’ ’ also 

to  measure  the  temperature-dependent  quenching  rate  coefficient  for  the 

v — R deactivation  process  in  DF(v.  = 1)  + HF(v  = 0)  collisions.  Because  of 

- 1 2 

the  large  endothermicity  (1055  cm  ) in  DF(v^  = 1)  + HF(v£  = 0)  collisions, 
the  v — v process  is  assumed  not  to  occur.  This  has  been  confirmed  in 
laser-excited  induced  fluorescence  experiments.  The  double  exponential 
decay  observed  in  HF(v^  = 1)  + DF(v^  = 0)  collisions,  for  example,  is  not 
observed  in  DF(Vj  = 1)  + HFIv^  = 0)  collisions.  The  upper  v^  levels  in 
DF(Vj)  + HF(v^  = 0)  collisions  also  reflect  v— ■ R transfer.  These  processes 
are  endothermic,  and  the  endothermicity  increases  with  increasing  v^  for 
single  quantum  v transitions.  We  assume  that  the  role  of  v— -R  energy 
transfer  in  HF(v^)  + HFfv^  = 0),  HF(v^)  + DFfv^  = 0),  and  DF(Vj)  + DF(v ^ = 0) 
collisions  can  be  deduced  experimentally  from  studying  v-*R  processes  in 

DF(v.)  + HF(v  = 0)  collisions.  Flow-tube  studies  have  been  carried  out  at 

* ^ 7 8 9 

room  temperature  by  Airey  and  Smith,  Kwok  and  Wilkins,  and  Smith7  on 

7 

v — R relaxation  of  DF(Vj  2 1)  + HF(v^  = 0).  Airey  and  Smith  reported  a 

36%  increase  in  the  v— ■ R rate  in  going  from  v.  = 1 to  v.  = 2.  Kwok  and 

8 2 1 1 9 

Wilkins  found  a v dependence  in  going  from  Vj=ltov^=3.  Smith  found 

a v dependance  in  going  from  v^  = 1 to  v^  = 5.  These  discrepancies  and 

lack  of  accurate  rate  data  for  upper  vibrational  levels  point  out  the  need  for 

further  experimental  measurements  and  theoretical  calculations  of  both  v 

and  temperature-dependent  rate  coefficients  for  the  energy-transfer  and 
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deactivation  processes  that  occur  from  upper  vibrational  levels  in  HF(Vj) 

+ DF(v^  = 0)  collisions  as  well  as  in  DF(v^)  + HFfv^  = 0)  collisions. 

Previous  trajectory  studieslO.il  on  HF(v^)  + HFtv^)  and  DF(Vj)  + 
DFfv^)  collisions  predicted  the  role  of  rotational  degrees  of  freedom  in 
reaction  dynamics.  In  particular,  it  was  found  that  one  or  more  quanta 
of  vibrational  energy  of  the  vibrationally  excited  HF  or  DF  molecule  is  con- 
verted into  rotational  energy  of  the  same  molecule  in  HF  + HF,  and  DF  + DF 
collisions.  This  conversion  of  vibrational  energy  into  rotational  energy  of 
the  same  molecule  occurs  by  an  intramolecular  energy-transfer  process 
and  is  responsible  for  formation  of  high  rotational  states  of  hydrogen  halides 
and  their  isotopic  analogs.  Such  high  rotational  states  cannot  be  formed 
from  chemical  pumping  reactions,  for  example,  such  as  H + F^  or  F + H^. 
The  v — v,  v-»R,  and  R — (R',T)  rate  coefficients  presented  in  these  pre- 
vious studies^'  H provide  the  rate  data  required  to  calculate  both  v and 
temperature-dependent  quenching  rate  coefficients  for  HF(Vj)  + HFfv^),  and 
DF(Vj)  + DFfv^)  systems.  The  v — R and  R — (R'.T)  rate  coefficients  were 
used  in  a rotational  nonequilibrium  model*  to  evaluate  the  temperature - 
dependent  quenching  rate  coefficients  for  self- relaxation  of  HF(Vj  = 1)  and 
DF(Vj  =1).  The  results  were  in  good  agreement  with  available  experi- 
mental data. 

A trajectory  study  was  undertaken  to  calculate  the  temperature- 
dependent  v — R rate  coefficients  that  occur  in  the  HF(Vj)  + DF(v^)  and 
DF(Vj)  + HFfv^)  systems.  This  work  completes  our  trajectory  studies  of 
HF(Vj)  + HF(v^)  collisions  and  their  isotopic  analogs.  The  v — R rate 
coefficients  are  used  to  calculate  temperature -dependent  quenching  v — R 
rate  coefficients  for  vibrational  relaxation  of  HF(v^  = 1)  by  DF(v^  = 0)  and 
DF(Vj  = 1)  by  HF(v^  = 0).  Rate  coefficients  are  calculated  for  12(v  — v) 
energy-transfer  processes  that  occur  in  HF(v^)  + DFfv^)  collisions.  Eleven 
of  the  twelve  v— v processes  have  not  been  experimentally  determined. 

jJj 

R.  L.  Wilkins  and  M.  A.  Kwok,  "Temperature  Dependence  of  HF(v^  = 1) 

+ HF(V2  = 0)  Vibrational  Relaxation,"  TR-0078(3603)-5,  The  Aerospace 

Corp.  , to  be  published. 
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H.  METHC"  OF  CALCULATION 


The  potential  energy  surface  used  here  for  HF-DF  and  DF-HF 
interactions  is  constructed  from  two  functions,  a London- Eyring-Polanyi- 
Sato  (LEPS)  potential  energy  function  for  the  short-range  interactions  and 
a point  charge,  dipole-dipole  potential  function  for  the  long-range  inter- 
actions between  the  four  atoms.  The  method  of  constructing  this  surface 
has  been  described**  previously  for  HF(v^)  + HFfv^)  interactions  and, 
therefore,  is  not  repeated  here. 

The  quasi-classical  procedure  described  in  Ref.  12  was  used  to 
examine  the  collision  dynamics  of  diatomic-diatomic  energy-transfer 
processes.  The  Hamiltonian  equation  that  describes  the  four-particle  sys- 
tem was  written  in  generalized  coordinates,  and  the  resulting  set  of  18 
equations  was  integrated  on  a CDC  7600  computer.  The  initial  vibrational 
states  of  HF(Vj,  J^)  and  DFfv^,  J were  assigned  the  vibrational  quantum 
numbers  v^  and  v^,  respectively.  Calculations  were  made  with  values  of 
Vj  and  v^  varying  from  0 through  4.  The  rotational  quantum  numbers 
and  J ^ were  Monte-Carlo  selected  for  the  calculations  involving  v — v and 
v — R energy-transfer  processes.  The  initial  values  of  the  coordinates  in 
each  collision  were  determined  by  means  of  (1)  a fixed  relative-translational 
energy  of  the  reactants,  (2)  a chosen  initial  relative  separation  between  the 
centers  of  mass  of  the  diatomic  molecules,  (3)  chosen  values  of  the  rota- 
tional quantum  numbers  and  J^,  (4)  the  impact  parameter,  (5)  the 
vibrational  phase  angles  of  both  molecules,  (6)  two  sets  of  orientational 
angles,  and  (7)  the  rotational  planes  of  both  molecules.  Averaging  over 

initial  internal  states  of  HF(v.,  J.)  and  DF(v_,J_)  was  carried  out  by  the 

11  L ^13 

technique  described  by  Porter,  Raff,  and  Miller.  For  excited  vibrational 
states,  the  calculation  of  the  initial  values  of  the  momenta  P^i  = 1 , . . . , 9) 
is  more  difficult  because  the  initial  amplitudes  of  both  molecules  are  no 
longer  the  classical  turning  points,  and  the  internal  momentum  vectors  are 
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no  longer  perpendicular  to  the  bond  axis.  For  this  case,  a method  based 
on  Euler  angles  and  adequately  described  in  Ref.  12  was  used.  The  distance 
R between  the  center  of  mass  of  the  HF  and  DF  molecules  was  taken  to  be 
8 a.  u.  This  value  is  large  enough  to  ensure  a negligible  initial  interaction 
energy.  Relative  translational  energies  were  assigned  values  that  ranged 
from  0.  5 to  6 kcal/mole.  Approximately  200  trajectories  were  run  for  each 
initial  set  of  parameters.  A uniform  distribution  was  used  for  the  square 
of  the  impact  parameter  b.  All  runs  were  made  with  a maximum  impact 
parameter  of  2.  5 A.  The  step  size  was  5.  7 X 10  *7  sec.  The  integration 
technique,  which  has  not  been  reported  in  the  literature,*  was  faster  than 
the  Runge-Kutta-Gill  procedure,  * ^ with  the  accuracy  of  the  integration 
tested  by  changes  in  the  step  size  and  by  integration  backwards  along 
selected  trajectories.  As  an  additional  verification  of  numerical  accuracy, 
all  trajectories  were  checked  at  each  point  along  the  trajectory  for  con- 
servation of  total  energy  and  total  angular  momentum.  Computation  time 
for  a single  trajectory  was  dependent  on  the  initial  parameters  but  was,  on 
the  average,  about  9 sec  per  trajectory  on  the  CDC  7600  computer. 

The  final  properties  of  each  trajectory  were  analyzed  to  determine 
the  nature  of  the  collision,  i.  e.  , the  total  angular  momentum  and  vibrational- 
rotational  energy  of  the  HF  and  DF  molecules.  The  partition  of  the  internal 
energy  between  vibrational  and  rotational  energy  was  determined  from  the 
internal  energy  of  the  HF  and  DF  molecule  and  their  total  angular  momentum. 
The  actual  technique  for  calculating  the  partitioning  of  the  internal  energy 
has  been  described  previously*^  and  is  not  repeated  here.  The  v— v energy- 
transfer  cross  sections  were  calculated  by  Method  2 described  in  a previous 
paper  by  Wilkins.  In  Method  2,  it  is  assumed  that  vibrational  energy  is 
transferred  in  every  collision  and  that  only  one  quantum  state  is  accessible 
in  the  transfer  process.  The  energy-transfer  cross  sections  and  specific 
rate  coefficients  were  calculated  by  means  of  the  equations  given  by 
Karplus,  Porter,  and  Sharma.  17 


G.  Emanuel,  personal  communication. 
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in.  RESULTS  AND  DISCUSSION 


A.  VIBRATIONAL- TO- ROTATIONAL  ENERGY  TRANSFER 

This  trajectory  calculation  indicates  that  collision  complexes  are  not 
formed  for  typical  HF-DF  or  DF-HF  collisions.  This  result  is  not  unex- 
pected inasmuch  as  the  potential  energy  surface  used  is  identical  to  that 
used  for  HF-HF  and  DF-DF  collisions.  Previous  trajectory  studies^'  H 
on  the  HF-HF  and  DF-DF  systems  indicated  that  dimer  formation  is  not 
required  to  explain  vibrational  self- relaxation  of  HF(Vj  = 1)  or  of  DF(Vj  = 1) 
in  the  temperature  range  at  and  beyond  300  K.  This  study  shows  also  that 
v — R processes  in  HF-DF  and  DF-HF  collisions  occur  by  intramolecular 
energy  exchange.  The  vibrationally  excited  HF  or  DF  molecule  can  convert 
one  or  more  quanta  of  its  vibrational  energy  into  rotational  energy.  This 
v — R mechanism  is  responsible  for  the  high  rotational  states  reported  by 
several  authors.  ^ Krogh  and  Pimentel‘S  results  indicate  that  the  direct 
population  of  high  rotational  states  by  chemical  pumping  reactions  can  be 
ruled  out.  The  v — R mechanism  proposed  by  Krogh  and  Pimentel^*  were 
observed  in  this  and  previous  trajectory  studies.^’  ^The  trajectory 
studies  predicted  rapid  v — R deactivations  with  large  AJ  changes.  Specific 
v — R rate  coefficients  for  HF-DF  collisions  are  given  in  Tables  1 through  3, 
and  for  DF-HF  collisions  in  Tables  4 through  6.  The  qualitative  results 
of  these  v — R processes  are  as  follows: 

1.  The  incident  vibrationally  excited  HF  or  DF  molecule 
converts  one  or  more  quanta  of  its  vibrational  energy  into 
rotational  energy  through  an  intramolecular  v — R energy- 
transfer  mechanism.  This  v — R mechanism  is  responsible 
for  the  formation  of  high  rotational  states  of  HF  or  DF. 

2.  The  energy  defects  across  the  reactions  for  these  v — R 
processes  have  smaller  values  than  would  have  been 
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Table  1.  Vibrational -to -Rotational  Detailed  Rate  Coefficients  for 
HF(v . = 1,  J = 2)  + DF(v  = 0)  at  T = 300  K 


Table  2.  Vibrational -to -Rotational  Detailed  Rate  Coefficients  for 
HF(Vj  = 2,  J = 2)  + DF(v  = 0)  Collisions  at  T = 300  K 
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(NJ 

(NJ 

(NJ 

(NJ 

(NJ 

(NJ 

(NJ 

(NJ 

(NJ 

(NJ 

(NJ 

(NJ 

(NJ 

(NJ 

(NJ 

(NJ 

(NJ 

(NJ 

(NJ 

(NJ 

(NJ 

(NJ 

(NJ 

(NJ 

(NJ 

(NJ 

in  units  of  lO1^  cm’/(mole-sec).  To  obtain  rate  coefficient  for  v -*  R relaxa- 
tion of  HF(Vj  = 2,  Jj  = 2)  by  DF(v£  = 0,  J^),  multiply  by  ( 2J ^ + 1) 

exp[-E(v2  = 0,  J2)/kT]/Q^(v2,  where  *8  t^ie  rotati°nal  partition 

function. 


Table  3.  Vibrational -to -Rotational  Detailed  Rate  Coefficients  for 
HF(v . = 3,  J = 2)  + DF(v  = 0;  Collisions  at  T = 300  K 
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expf-Efv^  = 0,  J^J/kTj/Q^Cv^,  J2).  where  Qr(v2>  * 2)  *8  l^e  rotati°na^  partition 
function. 


Table  4.  Vibrational -to -Rotational  Detailed  Rate  Coefficients  for 
DF(v,  = i,  J.  = 3)  + HF(v,  = 0)  at  T = 300  K 


(vl»  J v2#  vJf  Jj*  v2*  ^2^  kX  *s  tke  tota^  rate  coefficient  for  v -*  R relaxation 
in  units  of  10^  cm^ /(mole- sec).  To  obtain  rate  coefficient  for  v -*  R relaxa- 
tion of  DF(Vj  = 1,  J^)  by  HF(v2  = 0,  J^),  multiply  by  (2J^  + 1) 
exp[-E(v2  = 0,J2)/kT]/QR(v2,J2),  where  QR(v2>  J2  ) is  the  rotational  partition 
function. 


Table  6.  Vibrational -to -Rotational  Detailed  Rate  Coefficients  for 
DF(v . = 3,  J.  = 3)  + HF(v,  = 0)  Collisions  at  T = 300  K 


(vl’^l’v2’'l’^l,v2’^2^'  kX  **  rate  coefficient  for  v - R relaxation 

in  units  of  1012  cm3/(mole-  sec).  To  obtain  rate  coefficient  for  v - R relaxa- 
tion of  DF(v  ( = 3,  J^)  by  HF(v^  = 0,  J^),  multiply  k^  by  (2J^  + 1 1 
expl-Eiv^  = 0.  J2)/kT|/QR(v2,  J2),  where  QR(v2>  J2)  is  the  rotational  partition 
function. 


predicted  if  both  reactant  and  product  species  were 
assumed  to  be  rotationless. 

3.  Little  change  occurs  in  the  internal  energy  state  of  the 
target  molecule. 

Higher  rotational  states  of  DF  are  populated  in  DF(v)  + HF  collisions 
than  in  HF(v)  + DF  collisions  because  the  rotational  energy  of  DF  at  each 
value  of  J is  slightly  higher  than  one-half  that  of  HF  at  the  same  value  of  J. 
Since  a quantum  of  vibrational  energy  of  DF  is  approximately  75%  that  of 
HF,  the  J-states  of  DF  are  less  than  twice  those  of  HF.  The  net  effect 
indicates  only  a slight  increase  in  the  J-states  for  DF  as  compared  with 
those  for  HF.  The  results  in  Tables  1 through  6 indicate  that  the  deuterated 
species  of  the  hydrogen  halides  relax  more  slowly  than  their  hydrogenated 
analogs. 

Detailed  rate  coefficients  k , it,  / are  listed  in  Tables  7 

V1  v2  V1  v2 

through  10  for  v — R deactivation  processes  that  occur  in  HF(v^)  + DF(v^)  and 

DF(vj)  + HF(v^)  collisions.  The  overall  rate  coefficients  k^R  in  each  table 

increase  slightly  with  increasing  v.  Rate  coefficients  are  given  in 

Tables  8 and  10  for  collisions  of  vibrationally  excited  species,  i.  e.  , 

HF*  + DF*  or  DF*  + HF*).  Within  calculated  error,  k , ; /,  / and 

V1  v2  V1  v2 

k are  almost  identical.  Analysis  of  the  trajectory  results  indicates  a 
vK  0.  5 

temperature  dependence  of  T ' for  the  v — v and  v— ■ R energy-transfer 
processes. 

The  rate  coefficients  given  in  Tables  7 through  10  indicate  an 
increase  in  the  number  of  v — R channels.  Those  in  Tables  7 and  9 indi- 
cate a decrease  in  the  specific  rate  coefficients  for  v — R processes  with 
increasing  v for  collisions  involving  only  vibrationally  excited  species.  The 
overall  rate  coefficients  presented  in  Tables  7 and  9 are  not  comparable 
to  the  effective  quenching  rate  coefficients  measured  in  the  laboratory. 

The  predicted  large  increase  in  the  ratio  kvR(v  = 2)/kvR(v  = 1)  is  due  to  the 
increase  in  v-*R  energy-transfer  processes  that  results  from  the  opening 
of  two  channels  for  conversion  of  vibrational  energy  into  rotational  energy. 

g 

Kwok  and  Wilkins  found  that  the  rate  coefficients  for  v — R deactivation  of 
DF(v)  by  HF  have  a v^  vibrational  dependence.  The  rate  coefficients 
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Table  7.  Rate  Coefficients  kVj,  V2*  vf*  v?  for  Vibrational -to  - 
Rotational  Energy  Transfer  in  hE(v.  ) + DF(v  ) 
Collisions  at  T = 300  K 1 2 


V11 

v2 

i 

v2 

k , ,a 

vrv2;vrv2 

k a 
vR 

1 

0 

0 

0 

35.  7 ± 3.  0 

35.7  *3.0 

2 

0 

1 

0 

21.  1 ± 2. 7 

2 

0 

0 

0 

20.  2 ± 2. 7 

41. 3 ± 3.  1 

3 

0 

2 

0 

16.  7 ± 2.4 

3 

0 

i 

0 

15.  0 ± 2.  3 

3 

0 

0 

0 

15.4  ± 2.  3 

47.  1 ± 3.  1 

4 

0 

3 

0 

14. 5 ± 2. 3 

4 

0 

2 

0 

11.0  ± 2.  1 

4 

0 

1 

0 

12.  8 ± 2. 2 

4 

0 

0 

0 

11.4  ± 2.  1 

49.7  ± 3.  1 

i o o 

In  units  of  \.0VC  cmJ/(mole-sec). 
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Table  8.  Rate  Coefficients  kVj,  V2‘,  v[,  v£  for  Vibrational -to  - 
Rotational  Energy  Transfer  in  HF(v.)  + DF(v_) 
Collisions  at  T = 300  K 1 L 


V1 

v2 

vL 

k , / 

k a 

2 

Vv2;vl’v2 

VR 

1 

1 

1 

0 

1. 7 ± 0. 3 

1 

1 

0 

1 

2.  0 ± 0.  3 

1 

1 

0 

0 

1.  0 ± 0.  2 

4.  7 ± 0.3 

1 

2 

1 

1 

1.  5 ± 0. 2 

1 

2 

1 

0 

0.  9 ± 0.  2 

1 

2 

0 

2 

1.  8 ± 0. 3 

1 

2 

0 

1 

0.  5 ± 0.  1 

1 

2 

0 

0 

0.  7 ± 0.  2 

5.4  ± 0.4 

1 

3 

1 

2 

1.  3 ± 0.  2 

1 

3 

1 

1 

1.  0 ± 0.  2 

1 

3 

1 

0 

0.  6 ± 0.  2 

1 

3 

0 

3 

1.  5 ± 0.  2 

1 

3 

0 

2 

0.  7 ± 0.  2 

1 

3 

0 

1 

0.  3 ± 0.  1 

1 

3 

0 

0 

0.  5 ± 0.  1 

5.9*  0.  3 

2 

2 

2 

1 

1. 4 ± 0.  2 

2 

2 

2 

0 

0.  6 ± 0.  2 

2 

2 

1 

2 

1. 1 ± 0. 2 

2 

2 

1 

1 

0.  4 ± 0.  1 

2 

2 

1 

0 

0.  5 ± 0.  1 

2 

2 

0 

2 

1. 0 ± 0. 2 

2 

2 

0 

1 

1. 3 ± 0. 2 

2 

2 

0 

0 

0.  5 ± 0.  1 

5.  8 * 0.  3 

aIn  units  of  10^  cm^ /(mole-sec). 
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Table  9*  Rate  Coefficients  kvj(  y^;  v.»  v£  for  Vibrational -to - 
Rotational  Energy  Transfer  in  DF(v. ) + HF(v  ) 
Collisions  at  T = 300  K 1 2 


V1 

V2 

v2 

k a 

vi’v2;vl’v2 

k a 
vR 

1 

0 

0 

0 

32. 8 ± 3.  0 

32.8  ±3.0 

2 

0 

i 

0 

20. 2 ± 2.  6 

2 

0 

0 

0 

18. 5 ± 2.  6 

38.7  ± 3.  1 

3 

0 

2 

0 

14. 2 ± 2.  3 

3 

0 

1 

0 

10.9  ±1.9 

3 

0 

0 

0 

16.4  ± 2.4 

41.5  ± 2.9 

4 

0 

3 

0 

12.  3 ± 2.  2 

4 

0 

2 

0 

11.9  ±2.2 

4 

0 

1 

0 

11.  1 ± 2.  1 

4 

0 

0 

0 

12. 3 ± 2.  2 

47.  5 ± 3. 1 

aIn  units  of  10*2  cm^/fmole-sec). 
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Table  10.  Rate  Coefficients  kVj,  vj>  v£  for  Vibrational -to  - 
Rotational  Energy  Transfer  in  DF(v. ) + HF(vJ 
Collisions  at  T = 300  K 2 


V1 

v2 

v2 

k , /a 

vr  v2*  vr  v2 

k a 

VR 

1 

1 

1 

0 

1.  8 ± 0.  2 

1 

1 

0 

1 

1.9  ± 0.2 

1 

1 

0 

0 

0.9  ± 0.  2 

4.  6 ± 0.  3 

1 

2 

1 

1 

1.  6 ± 0.  2 

1 

2 

1 

0 

0.  9 ± 0.  2 

1 

2 

0 

2 

1.  6 ± 0.  2 

1 

2 

0 

1 

0.  6 ± 0.  2 

1 

2 

0 

0 

0.7  ± 0.  2 

5.4  ± 0.  3 

1 

3 

1 

2 

1.  5 ± 0. 2 

1 

3 

1 

1 

1. 0 ± 0. 2 

1 

3 

1 

0 

0.  7 ± 0. 2 

1 

3 

0 

3 

1.  5 ± 0. 2 

1 

3 

0 

2 

0.  6 ± 0.  2 

1 

3 

0 

1 

0.  2 ± 0.  1 

1 

3 

0 

0 

0.  5 ± 0.  1 

6.  0 ± 0.  3 

2 

2 

2 

1 

1.  5 ± 0. 2 

2 

2 

2 

0 

0.7  ± 0.  2 

2 

2 

1 

2 

1.  1 ± 0.2 

2 

2 

1 

1 

0.  4 ± 0.  1 

2 

2 

1 

0 

0.  5 ± 0.  2 

2 

2 

0 

2 

1. 0 ± 0.  2 

2 

2 

0 

1 

0.  3 ± 0.  1 

2 

2 

0 

0 

0.  5 ± 0.  2 

6.  0 ± 0.  3 

aIn  units  of  10^  cm^ /(mole-sec). 
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k / / for  deactivation  from  a given  initial  v.  level  to  the  final  v 

V1  v2-  V1  VZ  1 1 

level  are  approximately  equal.  The  results  in  Tables  7 and  9 imply  that 

the  total  quenching  rate  coefficient  for  v — R deactivation  should  perhaps  be 

partitioned  equally  among  the  (v^  — vj)  transitions  possible  from  a given  v 

level.  The  measured  v — R rate  coefficients  represent  total  deactivation 

rates  from  each  initial  vibrational  level  v.  The  product  hydrogen  halide 

molecules  are  distributed  over  all  the  vibrational  levels  below  the  initial 

vibrational  level. 

B.  VIBRATIONAL-TO- VIBRATIONAL  ENERGY  TRANSFER 

Because  of  the  large  endothermicity  (1055  cm  *)  in  DF(Vj  = 1)  + 

HF(v^  = 0)  collisions,  it  has  been  confirmed  in  laser-excited  induced 
fluorescence  experiments  that  the  v— -v  process  probably  does  not  occur. 

The  double  exponential  decay  observed  in  HF(v^  = 1)  + DF(v^  = 0)  collisions 
does  not  occur  in  DF(vj  = 1)  + HF(v^  = 0)  collisions.  Rate  coefficients  are 
given  in  Tables  11  and  12  for  room-temperature  v — »v  energy -transfer 
processes  obtained  from  this  trajectory  study.  The  collision  dynamics 
indicate  quite  clearly  that,  in  v— -v  transfer  in  HF(v^)  + DFfv^)  and  DF(Vj)  + 
HFfv^)  collisions,  vibrationally  excited  HF  molecules  transfer  vibrational 
energy  to  the  DF  molecules.  This  result  concurs  with  the  laser-excited 
induced  fluorescence  experiments.  The  v — v process  DF(Vj  = 1)  + 

HF(v,  = 0)-DF(v!  = 0)  + HF(v'  = 1)  has  a very  small  cross  section  because 
of  its  large  endothermicity  (1055  cm  ).  As  a consequence,  trajectory 
studies  were  completed  only  for  v— v processes  that  were  exothermic.  It 
was  not  possible  to  calculate  the  AE  energy  mismatch  for  the  endothermic 
and  exothermic  processes  without  the  appropriate  rate  coefficient  in  the 
endothermic  direction  for  the  v — - v process;  thus,  «.  olumn  6 of  Tables  11 
and  12  is  left  blank.  The  energy  mismatch  calculated  under  the  assumption 
that  both  reagent  and  product  HF  or  DF  molecules  have  zero  rotational 
energy  is  given  in  column  7.  Previous  trajectory  studies1®'**  have  shown 
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Table  12.  Rate  Coefficients  for  Vibrational -to-  Vibrational  Energy  Transfer 


that,  in  order  to  calculate  a realistic  energy-defect  for  these  processes, 

one  must  consider  the  rotational  states  of  the  reagent  and  product  species. 

This  trajectory  study  indicates  also  that  multiple  v— v exchange  processes 

do  not  occur  with  any  great  probability.  If  such  processes  did  occur,  we 

would  expect  v — v exchange  from  DF(v  = 2)  to  HF(v  = 1).  No  collisions 

were  observed  in  which  the  DF  molecule  transferred  vibrational  energy  to 

the  HF  molecule.  The  rate  coefficients  for  v— v processes  HF(Vj)  + 

DF(v^  = 0)  -*-HF(Vj  = Vj  - 1)  + DF(v^  = 1)  with  v^  = 1 through  4,  respectively, 

increase  with  increasing  vibrational  quantum  number  v..  The  computed  rate 

coefficients  for  the  v — v processes  were  k(v^  = 1,  = 0;  v|  = 0,  = 1)  = 

1.  6 X 1012  cm3 /(mole -sec)  and  k(v.,  v = 0;  v'  = 0,  v - 1,  v ' = 1)  = 

( v , - 1 ) l c 1 1 2 

1 • 4 1 k(l,0;0,  1)  at  300  K.  There  is  an  increase  in  these  rate  coeffi- 

cients with  decreasing  energy-defect. 

For  HF(v.  = 1)  + DF(v?  = 0)  collisions,  the  calculated  v — v room- 

vv  12  ^ 

temperature  rate  coefficient  k = (1.  6 ± 0.  9)  X 10  cm  /(mole-sec)  can  be 
compared  with  the  reported  (v  — v + v-*R)  values  of  2.  5 X 1012  cm3/(mole- 
sec)  by  Ahl  and  Cool3  and  Lucht  and  Cool,4  1.4  X 1012  cm3/(mole-sec)  by 
Bott  and  Cohen,  2 and  1. 2 x 1012  cm3/(mole-sec)  by  Hinchen.  1 It  is  diffi- 
cult to  assess  the  actual  value  of  rate  coefficient  for  v— v processes  because 
the  v— -R  and  the  v—v  processes  are  coupled  in  the  experimental  measure- 
ment. The  measured  rate  coefficients,  therefore,  represent  the  sum  of  the 
contributions  from  the  v — v and  the  v — R energy-transfer  processes.  The 
rate  coefficients  for  the  v — v processes  HF(v.  = 1)  + DF(v,)-HF(v'  = 0) 

+ df(v2  = v2  + *)  with  v2  = 1 through  3 decrease  with  increasing  v^. 

These  rate  coefficients  decrease  with  increasing  energy  defect.  The  rate 
coefficients  for  v — v processes  DF(v^  = 1)  + HF(v2>  — DF(vj  = 2)  + 

HFfv^  = v2  - 1)  with  Vj  = 1 through  3 decrease  slightly  with  increasing  v2> 
Since  there  is  an  increase  in  energy  defect,  we  believe  that  these  v—v  rate 
coefficients  should  exhibit  an  increase  with  increasing  v2.  The  rate  coeffi- 
cients for  v— v exchange  for  HF(Vj  = 1)  + DF(v2  = 1)  collisions  and 

DF(v.  = 1)  + HF(v.  = 1)  should  be  identical.  The  trajectory  study  predicts 
1 ^ 12  3 

a value  of  (7.  9 i 1.  8)  x 10  cm  /(mole-sec)  for  v—v  exchange  for 
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HF(vj  = 1)  + DF(v2  = I)  collisions  and  a value  of  (10.  1 ± 2.  0)  X 1012  cra^/ 
(mole-sec)  for  DF(v^  = 1)  + HF(v2  = 1)  collisions.  Both  values  are  within 
the  statistical  error.  This  type  of  discrepancy  might  explain  why  the 
calculated  v-»v  rate  coefficients  for  DF(Vj  = 1)  + HF(v2)  collisions  were 
found  not  to  increase  with  increasing  energy  defect.  There  are  no  experi- 
mental results  available  with  which  to  compare  these  rate  coefficients  for 
v—  v energy-transfer  processes  from  the  upper  vibrational  levels.  These 
theoretical  results  provide  insights  into  the  mechanisms  of  v — v energy- 
transfer  processes  from  upper  vibrational  levels  in  hydrogen  fluoride  and 
deuterium  fluoride  systems. 

C.  VIBRATIONAL  RELAXATION  OF  HF(v1  = 1)  BY  DF(v;  = 0) 

AND  OF  DF(v^  = 1)  BY  HF  (v..  = 0) 

The  detailed  rate  coefficients  for  the  v— R energy-transfer  processes 
listed  in  Tables  1 and  4 together  with  the  R — (R\  T)  rate  coefficients  pre- 
sented in  recent  papers  by  Wilkins  were  used  in  a rotational  non- 

equilibrium model  to  calculate  temperature -dependent  quenching  rate 
coefficients  for  vibrational  relaxation  of  HF(v.  = 1)  by  DF(v,  = 0)  and  of 
DF(Vj  = 1)  by  HF(v2  = 0).  A description  of  this  rotational  nonequilibrium 
model  was  discussed  previously  by  Wilkins  and  Kwok*  and  will  not  be 
repeated  here. 

In  Figs.  1 and  2,  the  results  are  seen  to  be  in  good  agreement  with 
measured  relaxation  rates.  Blauer  et  al.  22  reported  that,  within  their 
experimental  scatter,  the  vibrational  relaxation  times  of  these  two  systems 
are  identical  in  the  temperature  range  between  1500  and  4000  K.  The 
calculated  results  obtained  from  the  rotational  nonequilibrium  model  tend 
to  approach  those  of  Blauer  et  al.  22  The  relaxation  times  for  DF(Vj  = 1) 
by  HF(v2  = 0)  lie  between  those  for  DF(v^  = 1)  by  DF(v2  = 0)  and  HF(Vj  = 1) 
by  HF(v2  = 0).  Pt  is  larger  for  DF(v^  = 1)  + DFfv^  = 0)  collisions  than  for 

♦ 

R.  L.  Wilkins  and  M.  A.  Kwok,  "Temperature  Dependence  of  HF(v*  = 1) 

+ HF(v2  = 0)  Vibrational  Relaxation,"  TR-0078(3603)-5,  The  Aerospace 
Corp.  , to  be  published. 
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HF(Vj  = 1)  + HF(v2  = 0)  collisions.  The  v— R transfer  in  HF(Vj  = 1)  + 
HF(v2  = 0)  collisions  occurs  more  rapidly  than  in  DF(v(  = 1)  + DF(v2  = 0) 
collisions.  The  relaxation  times  for  HF'fv^  = 1)  by  DF(v2  = 0)  lie  below 
those  of  HF(v.  = 1)  by  HF(v_  = 0).  The  values  of  the  ratios  ku_  „ / 

kDF-DF'  ^F-HF^HF-DF'  and  kHF-HF^HF-DF  obtained  from  the 
theoretical  model  are  in  good  agreement  with  the  available  experimental 

data.  The  references  for  the  available  experimental  data  are  indicated  in 
Figs.  1 and  2 and  in  recent  papers  by  Wilkins.10,1 1,17  The  theoretical 
results  show  that  a model  that  considers  the  formation  of  high  rotational 
states  by  v— R intramolecular  processes  is  capable  of  explaining  the 
observed  temperature-dependent  quenching  rate  coefficients  measured 
between  300  and  1500  K. 


-30- 


IV.  CONCLUSIONS 


This  trajectory  study  predicts  that,  in  HF(v^)  + DF(v2>  and  DF(v^)  + 
HF(v2)  collisions,  v-»-v  processes  occur  by  inter  molecular  energy  exchange 
with  single-quantum  v transitions.  The  vibrationally  excited  HF  species 
transfer  vibrational  energy  to  the  DF  species.  The  vibrationally  excited 
DF  species  do  not  transfer  vibrational  energy  to  the  HF  species.  The  rate 
coefficients  for  the  v — v processes  HF(v.)  + DF(v  = 0)-HF(v'  = v.  - 1)  + 
DF(v2  = 1)  with  Vj  = 1 through  4,  respectively,  increase  with  increasing 
vibrational  quantum  number  v. . The  rate  coefficients  for  the  v— v 
processes  HFfVj  = 1)  + DF(v2>  — HF(v|  = 0)  + DF(v'  = v£  + 1)  with  v2  = 1 
through  3 decrease  with  increasing  v^.  In  the  v — R processes,  which  are 
nonresonant  energy-transfer  processes,  several  quanta  of  vibrational  energy 
HF(Vj)  or  DF(Vj)  are  converted  into  rotational  energy  of  the  molecule. 

The  v — R processes  occur  by  intramolecular  energy-transfer  mechanisms 
with  multiple-quanta  v transitions.  This  study  predicts  that  v— v up-pumping 
coupled  with  v — R deactivation  with  large  AJ  changes  provide  the  mechanisms 
to  produce  high  J states  in  HF  and  DF  chemical  lasers.  This  study  predicts 
that  collision  complexes  are  not  formed  for  typical  HF  + DF  or  DF  + HF 
collisions  at  300  K and  beyond.  A rotational  nonequilibrium  model  is  used 
to  calculate  temperature-dependent  quenching  rate  coefficients  for  vibra- 
tional relaxation  of  HF(Vj  = i)  by  DF(v2  = 0)  and  of  DF(Vj  = i)  by  HF(v2  = 0). 
The  results  were  found  to  be  in  good  agreement  with  available  experimental 
data. 
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THE  IVAN  A.  GETTING  LABORATORIES 

The  Laboratory  Operationa  ol  The  Aeroepace  Corporation  ia  conducting 
experimental  and  theoretical  inveatigatione  neceaaary  for  the  evaluation  and 
application  of  scientific  advances  to  new  military  concepts  and  systems.  Ver- 
satility and  flexibility  have  been  developed  to  a high  degree  by  the  laboratory 
personnel  in  dealing  with  the  many  problems  encountered  in  the  nation's  rapidly 
developing  space  and  missile  systems.  Expertise  in  the  latest  scientific  devel- 
opments is  vital  to  the  accomplishment  of  tasks  related  to  these  problems.  The 
laboratories  that  contribute  to  this  research  are: 

Aerophvsics  Laboratory:  Launch  and  reentry  aerodynamics,  heat  trans- 
fer, reentry  physics,  chemical  kinetics,  structural  mechanics,  flight  dynamics, 
atmospheric  pollution,  and  high-power  gas  lasers. 

Chemistry  and  Physics  Laboratory:  Atmospheric  reactions  and  atmos- 
pheric  optics,  chemical  reactions  in  polluted  atmospheres,  chemical  reactions 
of  excited  species  in  rocket  plumes,  chemical  thermodynamics,  plasma  and 
laser-induced  reactions,  laser  chemistry,  propulsion  chemistry,  space  vacuum 
and  radiation  effects  on  materials,  lubrication  and  surface  phenomena,  photo- 
sensitive materials  and  aeneora.  high  precision  laser  ranging,  and  the  appli- 
cation of  physics  and  chemistry  to  problems  of  law  enforcement  and  biomedicine. 

Electronics  Research  Laboratory:  Electromagnetic  theory,  devices,  and 
propagation  phenomena,  including  plasma  electromagnetics;  quantum  electronics, 
lasers,  and  electro-optics;  communication  sciences,  applied  electronics,  semi- 
conducting, superconducting,  and  crystal  device  physics,  optical  and  acoustical 
imaging;  atmospheric  pollution;  millimeter  wave  and  far-infrared  technology. 

Materials  Sciences  Laboratory:  Development  of  new  materials;  metal 
matrix  composites  and  new  forms  of  carbon;  test  and  evaluation  of  graphite 
and  ceramics  in  reentry;  spacecraft  materials  and  electronic  components  in 
nuclear  weapons  environment;  application  of  fracture  mechanics  to  stress  cor- 
rosion and  fatigue-induced  fractures  in  structural  metals. 

Space  Sciences  Laboratory:  Atmospheric  and  ionospheric  physics,  radia- 
tion from  the  atmosphere,  density  and  composition  of  the  atmosphere,  aurorae 
and  airglow;  magnetospheric  physics,  cosmic  rays,  generation  and  propagation 
of  plasma  waves  in  the  magnetosphere;  solar  physics,  studies  of  solar  magnetic 
fields:  space  astronomy,  x-ray  astronomy;  the  effects  of  nuclear  explosions, 
magnetic  storms,  and  solar  activity  on  the  earth's  atmosphere,  ionosphere,  and 
magnetosphere;  the  effects  of  optical,  electromagnetic,  and  particulate  radia- 
tions in  space  on  space  systems . 
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